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Abstract 



We investigate a class of models for new physics which could produce a large 
difference in sin 2/3 between J/ipKs and the "pure penguin" mode — > 

(j)Ks- In such models, the dominant effect is through a Z-penguin and therefore a 
pattern of deviation in sin 2f3 as measured in B^ Xi^s'^ VcKs', J/ip^S and ip'Ks 
is predicted. If the preliminary data concerning the discrepancy between J/ipRs 
and (pKs proves correct in magnitude, discrepancies in these other modes would 
likely be observable. We also consider the effects such new physics could have on 
the Bg system and the isospin analysis of B ^ Kir. We compare this scenario with 
a scenario where contributions to sin 2/3 in various modes is produced by gluino 
loops and down squark mixing. 
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1 Introduction 



The Standard Model (SM) provides a consistent explanation of CP violation which has 
been observed in the Kaon system and, more recently at the BaBar and BELLE B- 
factories. Indeed, the B-factories provide a number of decays where CP violation should 
be evident and therefore will allow stringent testing of the SM mechanism of CP violation. 

The most well established B-factory result is the determination of sin 2/? from the 
decay B ipKs and related processes P-jSj: 

sin 2P^K = +0.734 ± 0.054 (1) 

The SM expectation derived from the Kaon sector, the rate oi b ^ u transitions and the 
rate of BB oscillation is that |3]: 

0.64 < sm2pf^t < 0.84 (95% C.L.) (2) 

and so there is excellent agreement with the above result. 

A firm prediction of the SM is that the 6 — s penguin would have the same weak 
phase to order as the b ccs transition driving B ipKs, see |S] for experimental 
tests of the SM back ground. It therefore follows that in the "pure penguin" decay 
B (pKs the time dependent CP violation would be the same as in 5 — ipKs- Early 
results from the B-factories, however, seem to contradict this. If we interpret the time 
dependence of B ^ ipKs as a measure of /3, P^Ks^ the experimental results are: 

sm2p^Ks = -0.18 ±0.51 ±0.07 (BaBar^) 

sm2p^Ks = -0.73 ± 0.64 ± 0.22 {BELLE^) (3) 
sin 2/?<^K3 = -0.38 ±0.41 {Average) 

The result Eq. (jS)) provide a suggestive contrast to [3^Ks i^i Eq. ([T]). If this 2.7-sigma dis- 
crepancy proves, with more data, to be a real difference between sin 2(3^Kg and sin 2(5^Ks 
then it would provide definitive evidence for new physics (NP). 

This discrepancy is only possible if the NP contribution to B ^ 4>Ks is comparable 
to the SM. This places considerable constraints on the nature of the NP. In jH] it was 
argued that a promising class of models that can explain the effect involve flavor changing 
couplings of the Z-boson to bs. Such couplings may arise either from the rf-quark mass 
matrix if these quarks are mixed with extra vector-like down quark or through flavor 
changing penguin graphs containing new particles, such as the minimal supersymmetric 
Standard Model (MSSM) with large mixing in the up squark sector. Regardless of how 
the sZb coupling arise, this mechanism leads to a number of definite predictions. In 
this paper, we will discuss some of the implications of non-standard sZb couplings. In 
particular the same amplitude which gives rise to the 4>Ks anomaly should also contribute 
to related hadronic decays B — > Kit; and indeed to the golden modes such as ip^s and 
XcKs- The pattern of contributions to these case may provide additional evidence for 
this mechanism, and in general for any kind of physics beyond the SM. 
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While NP effects in various b qiq2(l3 decays have been estimated e.g. ^ we in- 
vestigate here differences in the CP asymmetry induced by interference between mixing 
and decay among different final states with the same ffavor content. If NP breaks parity 
and CP, we obtain in general different answers from mesons with different J'-'^ quantum 
numbers. 

In the case of i? — >• {cc)Ks where (cc) = J /ip^ip' ,rjc,Xi early experimental analysis 
have already been performed as part of the effort to determine sin 2(3 in the context of 
the SM. Their weighted average, in fact, enters the number given in Eq. (|T|) as 

siTi2P j^Ks{Ks^.+.-) = 0.82 ±0.08 
sm.2l3^i,iKs{Ks^-K+-K-) = 0.69 ± 0.24 

sin 2/3^^Xs = 1-01 ±0.40 (BaBar^) (4) 

sin 2/3^^X5 = 0.59 ±0.32 

Clearly at the present time sin 2/5 is consistent between these modes but it will be 
important to pursue further experimental study sensitive to the level at which NP enters 
the b CCS decay amplitude as well as in the corresponding S^-decays. Such studies 
can be done at the B-factories and/or at future hadron colliders. It is important to keep 
in mind that these searches for NP are based on the prediction of the SM that sin 2(3 is 
consistent for all such modes and is independent of fits to the unitarity triangle of sin 2/5 
which are currently dominated by a theory error of 14% that exceeds the experimental 
error in sin 2/3 by a factor of two, see Eqs. ((1} and 

The outline of the paper is as follows. In Section El we work out constraints on NP 
from B — s> 4>Ks data. In Section IHl the effect of non-standard Z-couplings in hadronic 
2-body decays are calculated and constraints are discussed. In Section 0] we compare 
the predictions of the scenario with non-SM Z-penguins with other models of NP and 
in Section El we conclude. In Appendix^ we give the matrix elements of B ^ (pK and 
B — > {cc)K decays in terms of the effective low energy Hamiltonian. In Appendix FBI the 
initial conditions and the leading log running of the Electro Weak penguins are specified. 



2 The B (j)Ks decay amplitude with NP 

Let us now assume that B (j^Ks receives a NP amplitude which accounts for the 
discrepancy in sin2(3^Ks- Thus, if a is the magnitude of the SM contribution and b is 
the magnitude of the NP contribution, then if we use a phase convention where the SM 
contribution is real the amplitude may be written as: 



A{(j)Ks) =a + be'^P+^^ 



(5) 

^Also Belle uses J/'i, rjc, xi in their analysis, however, do not give the individual contributions. 
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where p(A) is the strong (weak) phase difference between the two contributions and 
A{X) is the amphtude for ^ X while A{X) denotes the amphtude for S° X. 
The factor ri{X) denotes the CP eigenvalue of the final state X, e.g. ri{(f)Ks)=-l. 
It follows then that the time dependent CP asymmetry may be written as: 

T(B^(t) — » X) — T(B^(t) X) 
^cpit) = r(^o(^)^x)+r(i?o(t)^X) = cos(AmBt) + Sx sin(AmBt) (6) 

where the sine and cosine coefficients are given by: 

\A{X)\' - \A{X)\' 



C 



X 



\A{XW + \AiX)\^ 



2Im{A*{X)A{X)e-^^^) . _ 



2 



Here, 2/5 is the phase from the BB mixing defined in the phase convention of Eq.© 

From the observed time dependent CP violation of i? ^ (pKs we may obtain some 
information about these amplitudes on a model independent basis. In the following we 
will discuss the analysis of the data in the scenarios that the strong phase difference 
p = 0. 

In the case where there is no strong phase, then A*ri{X) = A and therefore Cx = 0. 
We can therefore extract the amplitude from the CP averaged branching ratio Fq = 
+ and Sx using: 



A{X) = ai 



ToiX)d-tr]{X)Sx + a2Jl-S_ 



e-'f" (8) 



where ai^2 = =tl giving a 4-fold ambiguity. Then, given (3 we can determine from Eq. Q 
for each value of a a value of b/a and A by 6/a e*'^ = A/ a — 1. 
Data are consistent with no direct CP violation in 5 — (pKs decay: 

C^K^ = -0.80 ± 0.38 ± 0.12 {BaBar[6l) 

C^Ks = +0.56 ± 0.41 ± 0.16 (BELLE^l) (9) 
= -0.19 ± 0.30 (Average) 

In Fig. d we show the weak phase A as a function of the NP to SM ratio b/a for 
the case without a strong phase. The curves are the solutions to Eq. (jHJ for the central 
value and the ± 1 sigma range of S^Ksi= sm2P^Ks) from the data average in Eq. © 
and (3 = 23.9°, the central value of the fit Eq. (j2)). The measured branching ratio 
B{B^ = 8.4±1.6- 10"^ [H] agrees within errors with the SM value, e.g. To 

avoid a "fine tuning" configuration i.e. a near cancellation between NP and SM we take 



TQ{(f)Ks)/4: < a < 2^TQ{(j)Ks). We recall that for fixed input S", Fq each value of b/a 
gives a 4-fold solution for A within and vr. Some solutions in Fig.Qend for larger values 
of b/a because of the lower bound on a. The cut-out regions around 6/a ~ 1 increase in 
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b/a 

Figure 1: 4- fold solutions following from Eq. (jS)) in the b/a - \ plane for S^^Ks = 
—0.79, —0.38, +0.03, i.e. the central experimental value and the 1 a range for vanishing 
strong phase. The curves corresponding to the central value (-0.38) are solid while the 
curves corresponding to the 1-sigma range (-0.79, +0.03) are dashed. For details see text. 



size if the upper bound on a would decrease and vice versa. We see that ratios b/a are 
typically order one with order one weak phase A, similar to the findings of Ref. [TU] . 

3 Z-penguin effects in 2-body 6-decays 

The Lagrangian of the effective FCNC sZb couplings maybe written as 

= f-2^r^ (hl.SLZs, + bn^.snZ:,) Z^ + h.c. (10) 
4:71^ 2 COS Ow ^ / 

where Z^h (2'^^) denote the left (right) handed coupling strength. The sZ6-couplings are 
experimentally constrained as 

^\Z,b + Z,bSM\^ + \Z'J^ <0.08 (11) 

which updates jTH [T2] , where details can be found. The bound in Eq. (fTT|) is based on 
inclusive B Xse~^e~ decays at NNLO JHl and corresponds to an enhancement of 2 — 3 
over the SM value ^ 

ZsbSM = -VtlVtsSm '' OwC^o ^ -0.04 , Z',, sm - (12) 

^We comment on NP in BB mixing in Section |S| 

■^The experimental bounds on the branching ratios of i? — > {Xs^ K, K*)£'^£~ decays have gone down, 
but the theoretical value has decreased from NLO to NNLO, too. 
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The inclusion of right handed couphngs into the distributions of semileptonic b — > si'^i' 
decay is straightforward since different hehcities do not interfere in the hmit of a vanishing 
strange quark mass. 

Tree level Z-boson exchange and subsequent qq pair production induces NP contri- 
butions to hadronic -B-decays via b —>■ sqq. These are proportional to the qqZ coupling, 
which we assume to be SM like, i.e., the coupling Zy to vector and axial vector Za 
currents are given, respectively, as Zy = hq — '^Qq sin^ 9w and Za = —hq- Note that we 
use Zyq^^^q + ZAq'y^'j^q such that the sign of Za is opposite to the PDG definition [T^ . 
The values and ratio of different qqZ coupling are compiled in Table^for sin^ 9w = 0.23. 





cc 


ss 


uu — dd 


Zv 


+0.19 


-0.35 


+0.54 


Za 


-0.5 


+0.5 


-1 


Za/Zv 


-2.6 


-1.4 


-1.9 



Table 1: Values of vector, axial vector Z-couplings and their ratio for different qq pairs. 

As can be seen from the ratio Za/Zy in Table [T] the cc system has the biggest spread 
among final states with different quantum numbers (vector versus axial vector coupling) 
induced by the Z-exchange. In the next Sections 13. Il and f3. 21 we estimate the implications 
of large, in general CP violating sZ6-couplings for B —>■ {cc)K and Bg — > {cc)(f) decays. 
Further, the couplings of Z-penguins to / = 1 mesons such as vr" are very large. In Section 
13. 31 we investigate whether current data on i? ^ Ktt^ decays yield additional constraints 
on the sZ6-couplings. Note that in order to explain a deviation in CP asymmetry from 
the SM as large as hinted by current data the couplings Z^'fj need to be near their upper 
bound and with a large phase. 

3.1 5-decays into charmonium 

We discuss the implications of non-standard Z-penguins in the decays B MK into 
charmonium M = r/c, ^' "5 Xo5Xi?X2- The NP effect in the decays is split according 
to the CP properties of the final (cc) states. Whereas the vector mesons couple to a 
vector current ~ Zy, the pseudoscalar and axial vector mesons couple to axial vector 
current ~ Za- Hence, the Z-penguin effect in rjc and xi is bigger than in the vector 
mesons " by a factor of 2.6, as can be read off Table ^ In the presence of a CP 
violating phase in the Z-contribution, this leads to a difference in sin 2(3 measured among 
the golden modes. The generic size of this effect is arg(P/T)P/T where P denotes the 
penguin contribution including the NP, and T is the SM tree contribution in the b — > ccs 
amplitude. In the SM, arg(P/T) ~ A^, whereas it can be order one in the presence of 
NP. Since B{B (j)K)/B{B — > {cc)K) ^ 10~^ an order one NP effect in 6 — sss gives 
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up to ten percent contribution to the b ccs amplitude. Therefore, we expect here 
differences in sin 2/3 of this order. 

To quantify this general prediction, we neglect for simplicity direct CP violation and 
obtain for the size of the NP effect 

sm2(3MKs - sin 2/3 = sin(2/3 - SiTg{MKs)) - sin 2/3 ^ -&xg{MKs) cos 2/3 (13) 



where we expanded in small phases in Aj A and abbreviated wgA{MKs)IA{MKs) = 
a.Tg{MKs)- Note that ri{xo,2Ks) = +1 and -1 for J/'^,'^'' ,ric,Xi- Numerically, we 
obtain in the non-standard Z-scenario 

0.11 forM = ^,^'■■■ 
|arg(Mi^<,)|< (14) 
0.25 for M = rjc, Xi 

Here we varied the magnitudes and phases of Zgb and Z^f, independently while respecting 
the constraint given in Eq. (fTTj) . These numbers as well as numerical estimates of other 
quantities to be discussed later in the text (also for the case with Z^j, = 0) are summarized 
in Table 121 Explicit formulae of the matrix elements and details can be seen in Appendix 
IXlandlEl Since 

aigiVKs) ~ Im{Ci + Cf + Ci + Cf + + ) ^ Zy ^ ^ 

the NP correction from Z-penguins has the same sign for vector mesons V = J/"^, as 
for the axial ones A = Xi^Vc- This correlation holds over the whole {Zgb, Z'^j^ parameter 
space after performing leading log QCD corrections. This is illustrated in Fig. |21 where 
the NP phases arg(y4i^'5') versus arg(\^i^'5') are shown for Z-penguins (x, blue). Also 
shown is the correlation in another NP scenario (+, green), the MSSM with additional 
flavor and CP violation in singlet down squark mixing, to be discussed in Section 14.11 
We estimate the Z-exchange effect in the difference between decays into axial and vector 
coupling charmonia as 

I sin2/?AXs - s^n2l3vKs\ < 0.18 cos 2/? (16) 

Available data given in Eq. Q are not significant yet i.e. the difference of the er- 
ror weighted average of axial minus vector coupling final states reads as sin2/3^xs — 
sm2f3vKs = -0.05 ±0.26. 

Decays into xo,2 final states are factorization forbidden modes since by C conjugation 
and Lorentz invariance 

(Xo(0++)|c7mc|0) = 0, (x2(2++)|c7;.c|0) = (17) 

and their amplitude A{B^ —>■ xo,2K^) requires gluon exchange. At order as the color 
suppressed Z-penguins compete with a color enhanced SM contribution from tree level 
W exchange, thus they are doubly l/Nc suppressed with respect to the factorization 
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Figure 2: The NP correction for axial coupling mesons aig^AKs) as a function of the 
one for vector coupling mesons aig{VKs) induced by non-standard Z-penguins (x, blue) 
and in the MSSM with down squark mixing 5j^R23 (+5 green) discussed in Section I^!T1 
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allowed modes. Therefore, the effect of NP from Z-penguins should be least pronounced 
in S — i> Xo,2Ks and time dependent CP asymmetries should return in these modes the 
least polluted value of sin 2/5 in the cc system. However, since naive factorization is badly 
broken in these modes, it is difficult to be quantitative here. Still, a departure in sin 2/5 
in xo,2 from the fitted value given in Eq. Q indicates the presence of NP, as well as any 
discrepancy in the extracted value of sin 2/5 among different charmonia. 

3.2 Z-penguin effects in the 5s-system 

We investigate here implications on 2-body decays in the S^j-system such as Bs 
{cc)(j) and Bg — >■ 00. The coefficient of the sm{AmB^t) term in the time dependent CP 
asymmetry reads as (see Eqs. © and ((Tj) with changes from B^, to S^-mesons) 

Sm^ = -v{M<j)) sin (arg M12 - &Tg{MKs)) (18) 

where we neglected for simplicity the width difference between Bg and Bg mesons and as 
before direct CP violation, i.e. we used |^/^| = 1. In the SM the phase of the mixing 
amplitude argM^^^ = —2/5^ = —2X^1] is tiny. Hence, we expect Scc(j, and S'^^ of 0{X^) 
in the SM. The sZ6-penguins do have a twofold effect in the -B^-system, on the decay 
amplitude similarly to the i^^-system and on the B^Bg mixing. Hence, both AB = 1 
and AB = 2 terms in 5*4/^ receive contributions from the NP. Concerning the latter and 
employing the notations of 15], we find 

M,% = ^f'^f, BsJlms^r^B (zt + Z'gf + 2Z:,Z%x) (19) 
SvT'^ sm Ow ^ ' 

Here, X = AP^^ = —2.84 [T^ includes differences in the bag factors (taken from lattice 
[17 ) and perturbative QCD corrections of the matrix element between operators with 
different Dirac structure, 7^L(i?) x 7^L(_R) and ^^L{R) x ^^^R{L). We find 

^ M£ ^ Aa ( Z:l + Z':^ + 2Zt,Z%X 

a correction of up to |r^| < 0.5 w.r.t. the SM. 

The analysis of NP on the decay amplitude is analogous to the discussion in Section 
13.11 However because the final state in these modes consists of two mesons of spiuT^ 0, 
to isolate the CP eigenstate components of the final state one has to perform angular 
analysis in Bg — > M0 decays with M = J/\l/, xi^2, 0- Otherwise, there is a dilution 
in SmiP from admixture of CP odd and even contributions which diminishes a potential 
NP effect and introduces an additional uncertainty. None the less, this is an excellent 
null test of the SM since any large CP asymmetry, even in the data summed over po- 
larization, would indicate the presence of NP. If the particle recoiling against the is a 
scalar/pseudoscalar, then of course there is a single amplitude and angular analysis is 
not required. This would be the case for xo and rjc- 
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With argMi2 = arg(l + rz) we find 

0.42 for M = ^, ^' •• 
I^Af^l < 0.47 foTM = ri„xi (21) 
0.66 for M = 

The spread in SM<t, between different charmonia is up to 0.15. The corresponding num- 
bers for the case with the flipped hehcity Z-coupling switched off for all quantities 
discussed in this section can be seen in Table |21 



3.3 Isospin analysis in B ^ Ktt decays 

We start with an analysis in i? ^ Ktt^ decays. Isospin analysis relates the decay 
amplitudes of neutral and charged B mesons as 



A{B^ ir°7r°) = B-A, - A{B+ K+TT^) =B + A (22) 

where B = B1/2 denotes the A/ = piece and A = yli/2 — 2^3/2 is a A/ = 1 combination 
of If = 1/2 and If = 3/2. Z-penguins do violate isospin and hence will give a significant 
contribution to A. Is the amount required to explain the anomaly in B <pKs i-e- of size 
comparable to the AJ = SM QCD penguins allowed by current data ? Experimental 
findings give 

- 1^ {B{B^ KH^) + B{B^ KH^))t{B+) I ~ ^■^'^ ^ ^-^^ 

This can be written in terms of a relative strong (weak) phase p(A) and ratio with 
relation 

cos A COS p = , ' ; ' (24) 

The 1 a allowed cos A, |v4/i?| parameter space for vanishing strong phases is shown in 
Figure El For non-zero p the allowed range increases and the constraint disappears for 
cosp = 0. Therefore, large isospin breaking contributions AjB ~ 0(1) are currently not 
constrained by S — *^ Ki:^ data if the phase A ~ 0(1) is large, if, for example, the NP 
contribution comes with a non CKM large weak phase. 

Another constraint on non-standard Z-penguins can come from B ^ tttt, B ^ Ktt 
data and SU{3) symmetry. These Neubert-Rosner type bounds constrain the ratio 

^ 2[BiB+^n'K+)+BiB-^7r'K-)] 

B{B+ TT+RO) + B{B- TT-Ro) ^ ^ 



in the presence of isospin and CP breaking NP by 

R-'> 



1 - e3/2V(l«l + I cos7|)2 + (|6| + I sin7|)^ 



(26) 
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Figure 3: Constraints on isospin violation from B —>■ Kn^ data for vanishing strong 
phases. The bounds weaken for non-zero strong phases. The range between the curves 
is allowed at 1 ex. 
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Figure 4: The ratio as a function of 7. Shown are the data at 1 a (solid), the bounds 
in a NP scenario with enhanced and order one CP violating Z-penguins (dotted) and 
the SM bounds for central values (dashed). 



where a, h are CP even (odd) isospin violating contributions, see [20] for details. There is 
a corresponding upper bound on obtained by interchanging the sign in front of the 
63/2 term. As discussed in Sectional semileptonic rare h — > decays bound enhanced 

Z-penguins to be at most 2 to 3 of their value in the SM, i.e. \a + ib\ < (2 — 3)|a5M|- 
Hence, if the SM is allowed, then the upper bound on is even less constraining for 
Z-penguins. However, the lower bound excludes large values of |a|, \b\, as can be seen 
from Eq. (j26p . For our analysis we find = 1.31±0.15, 63/2 = 0.20±0.02 using recent 
data [in] and use asM = 0.64, bsM = 1201 for the central SM value. 

We plot the ratio -R~^as a function of 7 in Fig. HI Displayed are the 1 a band 
from data (solid), the NP bounds with \a\ = \b\ = 2\asM\ (dotted) and for comparison 
the SM bounds (see [201) using central values (dashed). We see that for such moderate 
enhancement of the isospin violating contributions the NP bounds are even weakened 
compared to the SM, hence irrelevant. Note that in order to accommodate sm2f]^xs = 
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—0.4, one needs order one CP violation, i.e. \a\ ~ We checked that our conclusions 
are stable under variation of the sizes of a,b independently as long as \a\, \b\ < C(10). 

4 Predictions of Z-penguins vs other NP models 

In Sections 13.11 and 13.21 we have discussed the impact of complex, non-standard Z- 
penguins on time dependent asymmetry measurements B {cc)K and Bg — > (cc)0 
decays. In this Section we discuss further predictions of this NP scenario and compare 
them with those of other models. 

Non-standard Z-penguin effects in B ^ <pKs have been first discussed in [8^. Un- 
like B —>■ {cc)K decays where the chromomagnetic dipole operator is color octet and 
suppressed, NP in 6 — > sss decays can enter in both the 4-Fermi and the gluon dipole 
operators Ogj, see Appendix 1X1 and iBl We stress that the NP in a scenario with non 
standard Z-penguins is confined to the 4-Fermi operators, i.e. dominantly in the elec- 
troweak penguin Og. We checked that the NP effects induced in the dipole operators via 
leading log RGE mixing are below few percent. This is good since there are substantial 
theoretical uncertainties related to the matrix element of Ogj. Therefore, if it turns 
out that there is NP in B ^ (pKs decay but not in 6 — > ccs - something that can be 
checked by comparing sin 2/3 from decays into different charmonia - the scenario with 
non-standard sZ6-couplings - and any other who does not have CP violation beyond the 
SM in the chromomagnetic dipole operator - can be excluded. 

The correlation of CP asymmetries in i? — > (pKs with decays into charmonium with 
axial A = xi,Vc and and vector coupling V = J/\E', \E'' is shown in Fig. El To be specific, 
we plot sin2P AKs ~ sin2PvKs as a function of sin2P^Ks ~ sin2PvKs iii the Z'-penguins 
scenario (blue). Also shown is the outcome in the MSSM with down squark mixing 
^RR23 discussed in Section l^iTl for two values of the chromomagnetic as matrix element, 
Ks = —0.045 (red) and Ks = —0.030 (green). Data given in Eqs. ^ and (jH) yield 

sin2P^Ks ~ sin2/5(j/^^^/)i^g = —1.19 ± 0.42 (Average) (27) 

which is also displayed. The correlations shown in Fig. El from the two beyond the SM 
models are very distinct and allow to distinguish between them with improved data. 

4.1 Impact of super symmetric g — d loops on 6 ^ ccs decays 

Let us now consider the implications of the MSSM with additional flavor and CP violation 
beyond the Yukawa couplings on 6 — > ccs transitions. In particular, let us suppose that 
there is only mixing between the singlet scalar partners of the down quarks of 2nd and 
3rd generation, to be abbreviated here as S^j^23- This term induces FCNC though 
gluino-down squark loops. In this scenario, 6^j^ 23 the only source of beyond the SM CP 
violation and so this particular framework is rather predictive. For previous studies of 
gluino mediated effects in B (pKs see ||21j,|22|. 
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Figure 5: Difference in mixing [sin 2/5axs — sin 2/?i/Xs] ^ function of 

[sin 2(3^Ks ~ 2/5yxs] non-SM Z-scenario (blue) and in the MSSM witli 

additional flavor violation induced by S^^23- The latter is shown for Kg = —0.045 (red) 
and Kg = —0.030 (green). Also displayed is the 1 a range from data given in Eg. ()27|1 . 



The model gives contributions to the flipped 4-Fermi operators O3 g and the dipole 
operators Oy^, Og^. We use the initial conditions at the weak scale given in in the 
mass insertion approximation. We evolve the coefficients Ci from the SM and C'^ from 
the NP separately to the nib scale with the leading log RGE, see e.g [TH]. We scan over 
the squark mass ISOGeV < rriq < ITeV, the gluino mass 0.2 < {m-g/mqf < 1.3 and 
I^i?ij23l < 1 with arbitrary phase and include the constraints from data on B{h —>■ sj) 
following ^H]; pi]- Numerically, we find 

,,,,,,, 0.04 forM = ^,*'■■■ 
|arg(Mir5)|< (28) 
0.20 for M = rjc, Xi 

Similar to Z-penguins, the NP effects in decays to axial mesons A = Xiilc are bigger 
than in the decays to vector mesons V = J/\l/, We find that the relation 

s.rg{AKs) - +5.2 aTg{VKs) (29) 

holds up to a few percent as can be seen in Fig. |2l For the difference in sin 2/3 we obtain 

I sin 2(3aks - sin 2(3vks I < 0.19 cos 2(3 (30) 

The correlation with B 4>Ks is shown in Fig. for two values of the matrix 
element of the gluon dipole operator, Kg = —0.045 (red) corresponding to the asymptotic 

■^At large tan/? the magnitude of S^j^ 23 constrained by Bs A'^M^ data |25j . 
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distribution amplitude and Kg = —0.030 (green), the result in a quark model (see 
Appendix^. This illustrates the sensitivity to hadronic physics in this NP scenario. The 
difference in sin 2(3 between axial and vector coupling mesons is limited if the difference 
between (j)Ks and the vector ones is large and negative. Current data indicate that 
sin 2/3^7^3 ~ sm.2(3vKs ^ 0.04. We checked that this upper bound holds up to a value of 
the gluon matrix element which is 50 percent bigger than our central value Kg = —0.045. 

With values possible even greater than 60 ps~^ for Am^ the CP asymmetries 
in the i^^-system induced by ^^^23 '^^'^ dominate over the SM contribution 15.1ps~^ < 
Am^sAf < 21.0ps~^ @ 95 % C.L. 26 . Therefore, if at all measurable, this framework 
allows for large non-SM effects in S(^cc)<t> and S,j)^. 





both 




both 


z's, = o 


\^ig{V Ks)\ 


0.11 


0.07 


0.05 


0.03 


\8.ig{AKs)\ 


0.25 


0.18 


0.12 


0.08 


sin 2(3 AKs — sin 2(3vks 


0.18 cos 2/5 


0.12 cos 2/3 


0.08 cos 2/3 


0.04 cos 2/3 


sin 2(5^Ks - sin 2(3vks 


-0.47 ...+0.19 


-0.33.. .+0.14 


-0.20.. .+0.14 


-0.12. ..+0.10 


\rz\ 


0.50 


0.16 


0.19 


0.07 


\Sv<j)\ 


0.42 


0.08 


0.16 


0.02 


\^A(f)\ 


0.47 


0.17 


0.21 


0.04 


\S(f,^\ 


0.66 


0.42 


0.36 


0.17 


\SA<j> — Sv(j,\ 


0.15 


0.10 


0.07 


0.05 



Table 2: Upper bounds on and Bg quantities defined in text in the pres- 
ence of non-standard Z-penguins. For sm2P^Ks ~ sm2PvKs we show the accessible 
range. The first two columns are with the current constraints on the sZb couplings 
given in Eq. pi|l . whereas in the last two columns we entertain a scenario where 

\J\Zsh + ZshshA"^ + l^sbP — \ZsbSM\-, i-e., where data on B{h s£^£~) show no devi- 
ation from the SM. In the first and third columns both helicity couplings are present 
and the second and fourth ones are obtained with = 0. 



5 Conclusions 

One of todays most precise information on CP violation in the quark sector i.e. sin 2/? 
given in Eq. is an average over several final states with the same flavor content. While 
this procedure returns in the SM to very good approximation the same CP-asymmetries, 
it is wrong in general, for example if NP breaks parity and CP. 

We suggest here to study differences in sin 2/3 among decays B {cc)Ks into char- 
monia \E', r/c? Xo,i,2 with branching ratios of order 10~^ to search for physics beyond 
the SM. An order one NP effect m B ^ <pKs decay, which is required to explain the 
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current anomaly Eq. (jS)) or a similar large departure from the SM, leads to up to (9(10)% 
in the B {cc)Ks amplitude. This is just at the level at which the data agree with the 
SM, Eqs. and (j2|). We stress that differences in sin 2/5 can signal NP independent of 
improvements in the present error of the CKM fit. One might extend this program and 
compare CP asymmetries in decays into final states with the same flavor other than cc. 

We have explicitly shown that a NP model with non-standard Z-penguins does in- 
deed split sin 2/3 among different charmonia, see Eq. ()16|) . We stress that there is a strong 
correlation in this scenario between non-SM effects in 6 ^ ccs and b sss decays, since 
both get contributions from the modified 4-Fermi operators. If the (pKs anomaly per- 
sists, but sin 2/3(j/^_^/)Xc; — sin2/3(^j y,^)^^ vanishes, it model independently indicates the 
presence of an enhanced chromomagnetic dipole operator which carries a non-CKM 
CP violating phase. This emphasizes that b — >■ ccs and b sss decays are comple- 
mentary when constraining and distinguishing NP. A comparison of sin 2/? differences 
in the non-SM Z-scenario with the ones in the MSSM with gluino mediated FCNCs 
shows this in Fig. El We discussed time dependent studies in the corresponding decays 
of i?s-mesons, i.e. Bg — >■ (jxp and Bg — ^ [cc)(j). They show besides a similar effect from a 
NP weak phase on the decay amplitude one from NP in BgBs mixing and large effects 
are possible. The coefficient Scscf, in the CP asymmetry of Bg (cc)0 decays can be up 
to (9(0.4) in the presence of Z-penguins. In Bg 00 decay they induce an asymmetry 
up to order one, i.e. large as in 5 — (pKs decay. Our findings are summarized in Table 
El where we also entertain a scenario where the b — >■ si~^i~ branching ratio is SM like, 
hence the bound on the Z-penguins gets tighter \J\Zsb + Zgi, smI"^ + l^sfeP < \ZsbSM\ 
(assuming zero errors) than it currently is as given in Eq. (jllj) . There remain sizable 
effects in the observables discussed in this work, but closer to potential theory SM back 
grounds. Information on Z^2 phases also come from the Forward-Backward asymmetry 
m B ^ K*i~^i~ |T21 and inclusive B Xsi~^i~ decays in the long term future. 

We studied the effects of non-standard sZ6-couplings in an isolated manner. If this 
generic structure arises from a larger model beyond the SM e.g. the MSSM with addi- 
tional flavor violation beyond CKM, the signatures can be rather diluted. We recall that 
order one mixing between the cl and the squarks can generate via chargino-higgsino 
loops left handed sZ6-couplings near the experimental bound, whereas the right handed 
couplings are suppressed by the small strange Yukawa We show for comparison 

all quantities also with flipped helicity switched off in Table |21 The reach of the Z- 
couplings is sensitive to right-handed currents in particular in the i^^-system, i.e. Arris- 

In our study of CP violating time dependent asymmetries we neglected NP in transi- 
tions between the first and third generation down quarks which could affect BB mixing. 
Still, in this case, differences in CP asymmetries among different charmonia manifest 
the presence of new phases in 6 — ccs transitions. Others places where to look for 
Z-penguins are A;, — > A0 j2Hl, A^ A(cc), B {cc)K* and B — > (pK* decays where po- 
larization observables probe the handedness of the NP couplings. There is experimental 
support for the possibility of large electro weak penguins in B ^ Kit decays , which, 
for example, could be induced by non-standard Z-penguins. 

Acknowledgments G.H. would like to thank the theory groups at SLAC and LBNL for 
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A Matrix element oi B ^ {cc)K and B (f)K 

The effective Hamiltonian is given as, see e.g. [TH] 



with 



eff 



"71" 



Oi = {cLalt^bLI3){sLI3'ytiCLa) 

O3 = {sLlt^bi) J2 



i=l... 10,77,89 

02 = 
04 = 



{cLl^ibL){sLl^,CL) 



(A-1) 



O5 = (sL'Jf.bL) J2 ^Rl^^R 



O7 = -{,SLl^,bL) QqlRl^lR 



Oe = (sLa'jf.bLis) J2 (iRPl'^QRa 
3 

Os = -{sLalt,bLp)YQqqR(}l^qRa 



O9 = T^isLlt.bhjY.Ql'iLl^lL 



o 



77 



;mbSL(Jf,ubRF^'' 



O 



10 



Oh 



[sLali^bip) ^ QqQLpl^qLc 



m,SLaT:,a,,bRpG^^'' (A-2) 



167r2 



where the sum is over the active quarks q = u,d, s, c, b, and Qg is their electrical charge 
in fractions of e and a,f3 are color indices. The operators O^ are obtained from flipping 
L ^ R in the Oj. 

The matrix element of decays into mesons with vector coupling such as -B — J/ipK 
and B — ^ ip'K is proportional to {Nc is the number of colors) 



T = Vt,V,: (Oi + O3 + O5 + + + + Oy + O9 



c 



N, 



c 



The one describing decays into mesons with axial coupling such as i? Xi^ and 
B rjcK can be obtained from Eq. ()A-3|) by changing the sign of 05^6,7,8 and O^ 2 3 4 9 10 • 
The matrix element of 5 ^ (pK is proportional to 

+ K2C2 + nsCli^ +a^Cl 



(A-4) 

where ^2,8 stem from 0{as) corrections to the corresponding matrix element and the 
effective coefficient of the chromomagnetic operator, Cgg^ is defined e.g. in^H]- For NP 
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scenarios with an enhanced fesglue-couphng the matrix element of Og^ is particularly 
important. Keeping only terms leading in the heavy quark limit, we obtain 

in agreement with the QCD factorization jSUj calculation by [HI]. In the second step 
we evaluated the averaged momentum square from the gluon propagator 1/ < q^> with 
= ml{l — x) and convoluted it with the longitudinal light cone distribution amplitude 
of the meson, $||(a;), which encodes the momentum distribution of the constituent 
quarks in the meson, see [32] for details. Here, x denotes the fraction of the momentum 
carried by its s-quark. Theoretical uncertainties in Kg are from $|| and more generally, 
from power corrections to the factorization, which also limit the accuracy of the total 
amplitude in Eq. ()A-4|) . For asymptotic form ^\\{x) = 6(1 — x)x, which is supported by 
light cone sum rule calculations [32j we obtain m^/ < >= 3, bigger than the value from 
a quark model calculation ml/ < q'^>c^ 2. For our numerical study we use Kg = —0.045 
corresponding to the asymptotic distribution amplitude and K2 = —0.011 — i0.012 

B LLog renormalization for Z-penguins 

Z-penguins induce the following contributions at the weak scale 



C3 {m-weak) 



9' z:, 



sb 



9' z'x 1 



47r2 VtbVtl 6 



Cf (m_.) = - sin'^H.) 

Cf (m_,) = l^^^sin^^H. (B-1) 

Note that non-SM coefficients are treated as lowest order in and aw coupling con- 
stants. RG evolution to the mf,-scale is done with the effective 12 dimensional LLog 
anomalous dimension matrix, which is given in PH], except for the 7ifc, i = 7, . . . , 10 
and k = 77, 8g entries. They correspond to the Og-mixing of electroweak penguins onto 
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the dipole operators in the basis given in Eq. ()A-2|) and can be deduced from Plj as 

(7 = glKiQ^^m 

__16 _5 __H^ 

^7,77 ~ 9 ^7,83 - g 78,77 ~ 3X ^8.89 ~ 54 

_ 232 _ 59 _ 1180 g _ 46 

79,77 ~ "g]~ 79,8g ~ ~54 7lO,77 ~ ~g]~ 7l0,8g ~ ~27 

The sectors Oi and evolve independently with the same anomalous dimension. 
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